K — > 7T7T matrix elements beyond the leading-order chiral expansion * 
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We propose an approach for calculating K — > nn decays to the next-to-leading order in chiral expansion. A 
detailed numerical study of this approach is being performed. 



1. Introduction 

Controlling the effects of final state interactions 
(FSI) is one of the main barriers towards a high- 
precision theoretical prediction for K — > 7T7t de- 
cays. This is particularly difficult for lattice QCD 
because of the analytic continuation to Euclidean 
space. The finite-volume techniques developed 
in Refs. [|J3] can exactly take into account the 
FSI effects, but their numerical implementation 
is very demanding. In the foreseeable future, the 
most practical approach, which is also reliable 
and systematically improvable, to the study of 
non-leptonic kaon decays remains the combina- 
tion of lattice QCD and chiral perturbation the- 
ory (xPT). 

Apart from a study of the CP-conserving, AI = 
3/2, K — > 7T7T decay in Ref. all the exist- 
ing lattice results for matrix elements of the form 
{itit\Q\K) Q are obtained by simulating matrix 
elements of the kind (7r|Q|iif), then using lowest- 
order xPT H to relate them to the desired matrix 
elements. In this procedure, the effects of higher- 
order chiral corrections due to FSI, which could 
be very large, are completely missing. 
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3 In Ref. p|, the decay amplitude is also obtained at the 
precision of leading-order chiral expansion. 



Here we present the status of an on-going 
project in which all the relevant matrix elements 
for e'/e and the AI = 1/2 rule are being calcu- 
lated by computing matrix elements of the kind 
(ttit\Q\K), at some "unphysical" kinematics, on 
the lattice, and then using the chiral expansion 
to next-to-leading order (NLO) to obtain them at 
the "physical" kinematics. In this talk, we focus^] 
on the study of the chiral behaviour of AI = 3/2, 
K — > rnr decay amplitudes associated with the 
operators 

Qi = (8 a d a )L(upUp- dpdp) L + (s a Ua)L(updp)L, 

q=u,d,s,c 

3 ^ 

Qs = ^{s a dp) L ^ e q(q0a)R, 
q—u,d,s^c 

where a, f3 are colour indices and e q is the electric 
charge of q. $1^2) l,r means V>i7 M (l T 75)^2- 

2. Choice of "unphysical" kinematics 

The matrix elements (7r + 7r°|Qi|if +) (i = 4, 7 
and 8) are computed in the unphysical kinematics 
such that K + and one of the pions are always 
at rest, while the other pion might carry non- 
zero spatial momentum. We denote these matrix 

4 For details of other aspects of this work, please refer to 



2 



elements by 

(7r + (pV)7r°(0)|Qi|-ft: + )unphy S , if tt° is at rest, 
(7r + (0)7r°(p 7r )|(5 i |X + ) un p hys , if 7T+ is at rest. 

Notice that the spatial momentum p n in the 
above matrix elements might be zero as well. In 
this case, both pions are at rest. The correla- 
tors used to extract these matrix elements are 
discussed in Ref. ||. 

A technical difficulty arising at this stage 
is that the final states |7r + (p!„.)7r (0)) and 
|7r + (0)7r°(p 7r )) are not purely 1 = 2, because 
Bose-Einstein statistics does not rule out the 
1 = 1 components. In order to eliminate these 
components, we take the symmetric combination 

unphys 

(^+(0)^ (^)|Q l |^ + ) unphyf 



3. Chiral expansion at NLO 



(2) 



In order to relate ( 7r+7T °\Qi\K + ) s ^ >hys *° their 
counter parts at physical kinematics, denoted as 
(7r + 7r°|Q i |i ; r + )phys, we resort to NLO chiral ex- 
pansion. Under the chiral SU(3)t ® 5J7(3)r 
transformation, Q4 is in the (27,1) representation 
and (57,8 are in the (8,8) representation. There- 
fore the leading term in the chiral expansion for 
Q 4 (Q 7 ,s) is of 0(p 2 ) {0(p )). Chirally expanding 
Qi (Q 7 \s) to 0{p A ) (0(p 2 )) requires the knowl- 
edge of counterterm operators in the relevant chi- 
ral representation at this order, as well as one- 
loop calculations in xPT using the leading-order 
operators. 

There is only one X-PT representative for op- 
erators in (27,1) ((8,8)) representation at 0(p 2 ) 
{0{p )). At 0{p A ), there are thirty four %PT rep- 
resentatives for (27,1) operators, labelled as 
(i = 1, 2, . . . , 34) in Ref. 0. For the purpose of 
this work, we only need { 2 7) , 0^ 7) , of 7) , 7 27) , 
( 22 7) and [ 2 7) Q The seven 0{p 2 ) chiral repre- 
sentatives for (8,8) operators can be found in Eq. 

5 Notice that C^ 7) and ofP are not included in Ref. ||. 

(27') 

also contributes to (n\Q\K). 



(9) of Ref. @. Here we label them as 0\ ' ' 
(j = l,2,...,7). 

The one-loop xPT calculations at physical 
kinematics, in infinite volume, are reported in 
Ref. @ for (27,1), and Ref. § for (8,8). In 
our work, the chiral logarithms at the chosen un- 
physical kinematics are much more complicated 
than in these references, because of the energy- 
momentum injection at the operator. The re- 
sulting amplitudes depend upon three kinemat- 
ical variables, Mk > M„ (kaon and pion masses 
in the simulation) and E w (energy of the pion 
carrying a non-zero spatial momentum in the 
simulation), while they only depend upon Mk 
and M v (physical kaon and pion masses) at the 
physical kinematics. Furthermore, our numeri- 
cal simulations are performed in a finite volume 
in the quenched approximation. This requires 
one-loop calculations in finite-volume quenched 
%PT (qxPT). So far, we have only finished the 
one-loop (unquenched) x?T calculations in infi- 
nite volume as an exercise. These calculations 
have been performed in a general way, such that 
the four-momenta of K + . ir + and 7r° are used as 
the kinematical variables, and energy-momentum 
conservation is only implemented at the 0(p 2 ) 
strong vertices. The desired kinematics is im- 
posed at the end of the calculations. 

The results for the NLO chiral expansion for 
the amplitudes of interest are 



M 



(27) 



-6V2 



aM 2 - aM 2 



lphys " Ik f% " I * 
+ (p 4 - (3 5 +4/3 7 + 2(3 22 )M K 
+ {-40 2 + 3/3 4 + & - 2(3 7 
- 2(3 22 + 16/3 2 4)M 2 K M 2 

+ (4ft - 4/3 4 - 2(3 7 - l6fe 4 )M* 

a 1 , .(27) 
+ -T~n X (° ne lo °P)phy S ' 



M 



(8,8) 



2^2 



x i 7 



PhyS fKf 2 

+ [-{S 2 + S 3 ) + 2{S 4 + S 5 )+4S 6 }M 2 K 
+ [(S 1 +S 2 )+4(S 4 + 5 5 ) + 2S e ]M 2 \ 



3 



-(ram loop)gJ, 



MS 

for the physical kinematics, and 

(27)\sym _ ~6\/2~ 



(3) 



(M (Z '>) 



unphy S j R j2 



a(E 7r M^ 



E^M K , M^M K 



2 2' 
4/3 2 M^ + ( 4 /3 4 + 2(3 7 )E W M* 

(A - /?5 + 0t)M*M k 
(Pa -p 5 +P 7 + 2f3 22 )E 7I M 2 M K 
(-4/3 2 + 8/3 24 )m2m 2 



'2 

(-2/36 



+ (-16/3 24 )^M 



4/3 7 + 4f3 22 )E^M n M K 
%)M K 
1 1 2f3 22 ElM v M K 



+ (ft + 2/3 7 ) (M„ + ^ W )M 3 



+ 



x [(one loop)^)]^ 



unphys' 



(^(8,8)) 



sym _ 



2V2 



x 1 7 



Unphys /k/ 2 

+ [4(* 4 + <5 5 ) + 2<5 6 ]M 2 + [_($ x + 5 2 )]^M W 
+ [~(*i + fc) - («a + 5 3 )](M 7r + E^)Mk 



[2(6 4 + S 5 ) + 4S 6 ]M 2 K 



fx ft 



[(one loop) 



(8,8)isym 



J unphys' 



(4) 



for the unphysical kinematics described in Sec. 
H Here a (7) is the coupling constant accom- 
panying the 0(p 2 ) (O(p )) operator, and ft (Si) 
are the coupling constants associated with 0(p A ) 
(0(p 2 )) counterterm operators for (27,1) ((8,8)) 
chiral representation. We are able to reproduce 
the (one loop)^; s and (one loop)p h ' yi ; results in 
Refs. jl0| and l^j. We have also checked that all 
the amplitudes in Eqs. (||) and (^) are indepen- 
dent of the renormalisation scale p x in x?T with 
the same p x dependence in each NLO coupling. 

By having enough data points for the ampli- 
tudes in Eq. (^) at different values of E n , and 
Mki one can determine the coupling constants 
up to NLO chiral expansion and construct the 
physical amplitudes in Eq. (0). We are currently 



performing a detailed numerical study (clover ac- 
tion, j3 — 6.0) of this strategy, as well as the one- 
loop calculations in finite- volume q%PT. Figure 1 
shows an example of our numerical work. 




(aM K ) 

Figure 1. Bare lattice (tttt(I = 2)\Qg\K), in lat- 
tice units, plotted against (aMx) 2 ■ Circles are 
the data points in which both pions are at rest. 
Diamonds are the data points in which one of the 
pions carries non-zero spatial momentum. Points 
with the same (aMK) 2 are shifted for clarity. The 
solid line is the matrix element in the chiral limit, 
obtained via a Rt to Eq. (Q) using only data 
points where all the energies are below 1 GeV 
and setting chiral logs to zero. The dashed lines 
are the associated statistical errors. 
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